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ABSTRACT 

Distributed parameter modeling is being seen to offer a viable alternative to finite element 
approach for modeling large flexible space structures. The introduction of the transfer matrix 
method into the continuum modeling process provides a very useful tool to facilitate the distributed 
parameter model applied to some more complex configurations. A uniform Timoshenko beam 
model for the estimation of the dynamic properties of beam-like structures has given comparable 
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results. But many aeronautical and aerospace structures are of the non-uniform sections or 
sectional properties, such as aircraft wing, satellite antenna. 

This paper proposes a piecewise continuous Timoshenko beam model which is used for the 
dynamic analysis of tapered beam-like structures. A tapered beam is divided into several segments 
of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite element 
analysis, the closed-form solution of the Timoshenko beam equation has been used. Application 
of transfer matrix method relates all the elements as a whole. By corresponding boundary 
conditions and compatible conditions a characteristic equation for the global tapered beam has been 
yielded, from which natural frequencies can be derived. A computer simulation is shown in this 
paper, and compared with the results obtained from the finite element analysis. While piecewise 
continuous Timoshenko beam model decreases the number of elements significantly, comparable 
results to the finite element method are obtained. 

SYMBOLS 

A sectional area, or characteristic matrix 

a parameter in the Timoshenko beam equation, a 2 = El/m 

C i ,C2,C3,C4 mode shape coefficients 
Det [A] characteristic determinant 

E modulus of elasticity 

G shear modulus 

I second moment of area of the beam section 

k bending stiffness, k = El 

Lj length of the jth beam segment 

M bending moment 

m mass per unit length of the beam 

Q shear force 

r radius of gyration of the beam section, r 2 = 1/A 

T time function 

t time 

x, y, z Cartesian coordinates 

y(z,t) lateral deflection 

z dimensionless z-coordinate, z = z/L 

Y (z) spatial lateral deflection function 

a, y dimensionless parameters 

P eigenvalue coefficient, P 4 = co 2 / a 2 
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£ Timoshenko shear coefficient 

V slope of lateral deflection 

to circular natural frequency 

[$] transfer matrix for the global beam 

[<I>]j transfer matrix for the jth beam segment 

<I> id elements of the global transfer matrix 

<)> id^ elements of the jth beam transfer matrix 


1. INTRODUCTION 

Distributed parameter modeling is being seen to offer a viable alternative to finite element 
approach for modeling large flexible space structures. Continuum models have been made of 
several flexible space structures, which include the Spacecraft Control Laboratory (SCOLE) [i], 
Solar Array Flight Experiment p], NASA Mini-Mast Truss pj, the Space Station Freedom [4]. 
Especially, the introduction of the transfer matrix method into the continuum modeling process 
provides a very useful tool to facilitate the distributed parameter model applied to some more 
complex configurations [5,6]. A uniform Timoshenko beam model for the estimation of the 
dynamic properties of beam-like structures has given comparable results [7]. But many 
aeronautical and aerospace structures are of the non-uniform sections or sectional properties, such 
as aircraft wing, satellite antenna. 

This paper proposes a piecewise continuous Timoshenko beam model which is used for the 
dynamic analysis of tapered beam-like structures. A tapered beam is divided into several segments 
of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite element 
analysis, the closed-form solution of the Timoshenko beam equation has been used. Application 
of transfer matrix method relates all the elements as a whole. By corresponding boundary 
conditions and compatible conditions a characteristic equation for the global tapered beam has been 
yielded, from which natural frequencies can be derived. A computer simulation is shown in this 
paper, and compared with the results obtained from the finite element analysis. While piecewise 
continuous Timoshenko beam model decreases the number of elements significantly, comparable 
results to the finite element method are obtained. 


2. TRANSFER MATRIX OF A TIMOSHENKO BEAM 


Timoshenko beam model accounts for both rotary inertia and shear deformation of the 
beam. Usually, Timoshenko beam model produces more accurate estimation of the modal natural 
frequencies compared with the Bemoulli-Euler beam equation, especially for the range of higher 
frequencies [8]. In this section, a transfer matrix for Timoshenko beam model has been derived. 
The Timoshenko beam is represented by the equation, 

jjL + xa.^-1 . _QL ( 1 +-E_) 3y _+ m 2 =o (2.1) 

az 4 m at 2 n* eG az 2 at 2 eega 2 at 4 


For harmonic motion, y(x,t) can be expressed as 


y (x,t) = Y(x) e J “* 

then, Eq.(2.1) will become 

(l + ^_)(o2Y'+(-m^- - £ )Q )2 y =o (2.2) 

EA eG eEGA 2 EI 

Defining P 4 =co 2 /a 2 , where a 2 =EI/m, Eq.(2.2) becomes 


Y" + p 4 r2 (1+ Jt)Y" + p 4 (pV(i)-l]Y =0 


(2.3) 


where, r 2 =I/A, the radius of gyration of the section. If we use the following dimensionless 
parameters. 


z = 


L’ 





the Timoshenko equation may finally be written as 


Y" + (PL) 4 (a + Y) Y* + (PL) 4 [ (PL) 4 ay- 1 ] Y = 0 (2.4) 


Assuming that the solution is 


Y(z) = Ae( BL >* 


which, when substituted into Eq.(2.4), leads to 



(2.5) 




i tu 


I L 



(BL) 4 + (PL) 4 ( a + y ) (BL) 2 + (pL) 4 [ (PL) 4 ay - 1 ] = 0 
The solution to the Eq.(2.5) is as follows. 

(BL),, 2 = ± OIL) = ±S [ . (PL)"(a + Y) + V(PL) 8 (a-Y) 2 +4(PL) 4 
(BL) 3 , 4 = ± j (BL) = ± j & [ (PL) 4 (a + Y) + V (pL) 8 (a - y) 2 + 4(PL) 4 ]‘' 2 

Then, the solution to the Eq.(2.4) can be expressed as 

Y(z) = Ci sin (0Lz) + C2 cos (0Lz) + C3 sinh (r|Lz) + C4 cosh (r|Lz) (2.6) 

Similarly, for the bending slope the differential equation has the same form as the 
Eq.(2.4), 

+ (PL) 4 ( a + y) + (pL) 4 [ (pL) 4 ay - 1 ] ¥ = 0 (2.7) 

The solution to Eq.(2.7) will be 

'P(z) = CT1C1 cos (0Lz) - CT1C2 sin (0Lz) + O2C3 cosh (r|Lz) + sinh (r|Lz) (2.8) 
where. 


(eL)--O-(pL) 4 

and a 2 = ;L 

^ L ) + 7^T^ L > 4 

(0L) J 

L 

L Cn L ) 


For the Timoshenko beam model, the shear force is 


3 3 

Q(z,t) = k-L -k 

8 z 3 eGA 3z3t 2 


JEL 


a 3 y 


at 2 


or, equivalently. 


And the bending moment is 


Q(z) = k Y"+ k ^ y'+ Jgj 2 y 
eGA 


M(z,t)=k^ -k-m- 
0 Z 2 eGA 



or, equivalently, 


(2.9) 


( 2 . 10 ) 

( 2 . 11 ) 
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( 2 . 12 ) 


M(z) = kY%k^ y 

eGA 

Eqs. (2.10) and (2.12) can be written in dimensionless format as, 

Q(z) = ^- Y"+ ^a(pL) 4 Y + ~t(PL) 4 V (2. 1 3) 

M(z) = -k- Y'+ J^X(pL) 4 Y (2.14) 

L 2 L 2 

Substituting the solutions of Y(z)(Eq.2.6) and 4 , (z)(Eq.2.8) into Eqs.(2.13) and (2.14), we derive 
Q(z)=-kanCicos0Lz+kanC2sin0Lz+ka2iC3Coshr|Lz+ka2iC4sinhriLz (2.15) 

M(z)=-kOi2Cisin0Lz-kai2C2Cos0Lz+ka22C3sinhriLz+ka22C4CoshriLz (2. 16) 


where, 

an = 0 3 - aL 2 p 4 0 - 7 L 2 p 4 ai ? 
a 2 i = r| 3 + aL 2 p 4 T] + 7L 2 p 4 a 2 ? 


an = 0 2 - aL 2 p 4 
c 2 2 = r\ 2 + aL 2 p 4 


For the jth beam element, the displacement Y(0), slope 'F(O), shear Q(0), and the bending 
moment M(0) at the end of z=0 can then be written in matrix form as, 


Thus, 


where. 


’ Y 0 1 








' 0 

1 

0 

1 ' 

(C i) 



<*i 

0 


0 

)c 2 

Qo 


■ka n 

0 

kCT 2 l 

0 

H 


0 

-kai2 

0 

ka22. 

j lc 4 /j 

, Mo 

j 







Y 0 

Qo 

Mo 




- j 


(2.17) 


(2.18) 
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L /"k 


O l 0 CT2 0 
-koii 0 k02i 0 

0 -k<Jl2 0 k022 Jj 


012+022 


012+022 


O1O21+O2O11 


Oil 

01021+O2OII 


k(0i021+020n) 


Ol 

k(0i02i+020n) 


k(0i2+022) 


k(0i2+022) 


Similarly, at the end of z=L, the corresponding quantities are, if written in matrix form. 


where, 


T La, Is! 


ic 4 i 


Mj- 


sinQL cosQL sinhrjL coshriL 

aicosGL -ctisinGL a 2 Coshr|L a 2 sinhr|L 

-kcncosSL ko n sin0L k<72icoshr|L ka 2 isinhT|L 

-kOi2sin0L -kGncosGL ka22sinhr|L ka22Coshr|L Jj 


Substituting Eq.(2.18) into Eq.(2.19) we obtain 


* L Uwj r ° 

Ql Qo 


M l J j 


Mo I j 


where, [<t>]j is the transfer matrix of the Timoshenko beam. 


and the elements of the transfer matrix [d>]j are as follows. 


CPl I 1 (CT22COS0L + ai2C0ShT|L) 

012+022 


(pi 2 = 1 (G2isin0L + CTi isinhr|L) 

Cia2i+cr 2 an 
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<Pl3 

9l4 

921 

922 

923 

924 


X 


k(aia 2 i+a 2 On) 


(- c 2 sin0L + Cisinh'nL) 


X 


k(ai 2 +a 22 ) 


(- cosGL + coshriL) 


X 


(aiO 2 2Sin0L - CT2<7i 2 sinhr|L) 


Oi 2 +a 2 2 

^ ((Ti<J 2 icos0L + a 2 ancoshT|L) 

aia2i+(J2<7u 

gi c 2 


k(<3ia 2 i+a 2 an) 


(- cos0L + coshr|L) 


X 


k(oi 2 +o 22 ) 


(<7isin0L + <r 2 sinhTiL) 


(p 31 = k (aiia 22 sin0L + <Ti 2 a 2 isinhriL) 

ai 2 +a 22 

<p 32 = — ( C0s e L . coshtiL) 
<Tia 2 i+a2<Tn 

q, 33 = 1 (o 2 ancos0L + aia 2 icoshTiL) 

<7i<T 2 i+<j 2 an 

934 = 1 (- <Tysin0L + a 2 isinhr|L) 

Ci 2 +<T 22 


CP 41 = - (cos0L - coshriL) 

G12+G22 

942 = - — - j — (ai 2 a 21 sin0L - anC 22 sinhTiL) 

®1<J21 + ®2<*11 

943 = — — f~ - (a 2 o i 2 sin0L + aia 22 sinhT|L) 

aia 2 i+a 2 an 

944 = ^ (ai 2 cos0L + a 22 coshr|L) 

Oi2+0 22 


3. PIECEWISE CONTINUOUS MODEL FOR A TAPERED BEAM 

A tapered beam can be considered as a piecewise continuous step beam consisting of N 
uniform beam elements as shown in the figure. Using the transform matrix (Eq.2.21) to describe 



each beam element, then the state vectors at the two ends of the global beam will be related by the 
global transfer matrix [<&], that is, 
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(3.1) 




E=:3 


S3 



Y 


V 


-Mi 

¥ 

Q 

Q 

M 

N 

M 


where, the global transfer matrix 

w=riWj 

j=N 


Without loss of generality, let us consider the case of N=3. As N=3, the global transfer matrix 
will be 



(3.2) 


where, the superscripts (j) represent the jth beam, and <j>ia(i) is the elements of the transfer matrix 
for the jth beam. 

If we consider a cantilevered beam fixed at the end of z=0, we have the following boundary 
conditions: at the fixed end: Y(0)=0 and \|f(0)=0; at the free end: Q(L)=0 and M(L)=0. Applying 
the boundary conditions to the global equation (3.1), we will have 


Y 

to\ 

* 

=[*] 0 

0 

Q 


\ M ) 

1 0 

3 


(3.3) 


Rearranging the state vector, Eq.(3.3) can be written as 


44 


= [ 0 ] 


(3.4) 


where, the coefficient matrix 


- 1 0 < 1*1 3 4*14 

0 - 1 4>23 4*24 

0 0 4*33 4*34 

0 0 4*43 4*44 


and <3>jj's are the elements of the global transfer matrix. The condition for Eq.(3.4) to have a non- 
trivial solution is that the determinant of the coefficient matrix equals zero, that is 

- 1 0 <Di3 4»m 

Det[A] = Det 0 ‘ 1 4,23 4,24 =0 (3.5) 

0 0 4*33 4>34 

_ 0 0 4*43 4*44 - 

Eq.(3.5) is the so-called characteristic equation. Solving for the roots of the characteristic 
equation, we obtain the natural frequencies CD's. Expanding the determinant in Eq.(3.5) we can 
simplify the characteristic equation as 

<1*33 <1*44 * <1*34 <1*43 = 0 (3-6) 

or, expressing Eq.(3.6) in terms of the elements of each sub-transfer matrices, we have 

&(£ **) 4 # < 3 - 7> 

If we consider a free-free beam, then the boundary conditions will become as Qo=Mo=0 and 
Q3=M3=0 at the both ends. Thus the characteristic equation (Eq.3.5) should be 


Det [A] = Det 


4*31 4*32 0 0 

4*41 4*42 0 0 

4>ii 4>i2 -1 0 

4*21 4*22 0 -1 


45 


<I>31 ^*42 ' ^*32 ^*41 = 0 


( 3 . 9 ) 


Expressing Eq.(3.9) in terms of the elements of each sub-transfer matrices, we have 


I (i [| (t At?j*8'] ■ [i (| [| (| Af?)*!! 1 ] = 0(3.10) 


4. COMPUTER SIMULATION 

The computer simulation is designed to analyze the natural frequency for a tapered beam 
with 15-meter length (Fig.4.1). The modulus of elasticity is assumed to be E=200*10 9 N/m 2 . To 



simplify the calculation, the change of the sectional foil is specified by the changes of the second 
moment of area of the beam section and the mass of the beam segments along the longitudinal axis 
z, that is, assuming 

I = (0.02z 2 -0.6z+5.375)* 10' 6 (m 4 ) (4.1) 

and 

m = 0.0 1 1 2z 2 -0.495 z+7.708 (kg.) (4.2) 

In so doing, we may readily determine the sectional parameters necessary for the element 
stiffness and mass matrices when we divide the beam as any desired number of segments. For 
example, if we use three uniform beam elements to represent the global tapered beam, then we use 
Z!=2.5, Z2=7.5 and z 3 =12.5 to calculate the Ij and mi for each beam element according to Eqs.(4.1) 
and (4.2). They are 



Ii=4*10' 6 , I 2 =2*10‘ 6 , I 3 =l*10- 6 (m 4 )andm 1 =6.54, m 2 =4.03, m 3 =3.27 (kg) 

In the computer simulation, three-segment piecewise continuous Timosenko beam model 
has been applied. For finite element analysis, the commonly used two-node and four-degree-of- 
freedom plane beam element has been selected. Table 4.1 exhibit the comparison of the frequency 
results calculated by both the piecewise continuous Timosenko beam model and the finite element 
model. The results show that at least ten beam elements are needed for the finite element analysis 
to achieve the comparable frequency values while the piecewise continuous Timosenko beam 
model uses only three beam segments. The advantage is clear in decreasing the number of 
elements by using the piecewise continuous Timosenko beam model to analyze large flexible 
tapered beam-like structures. 


Table 4.1 The comparison of the results obtained from 
Finite Element Method & Piecewise Continuous Timoshenko Beam Model 

( Circular Natural Frequency, rad/sec ) 


Order of Mode 

1 

2 

3 

4 

5 



3 

15.274 

70.066 

179.455 

415.646 

775.990 

.9 

X 

4 

13.364 

61.444 

159.541 

305.038 

577.351 


§ 

a 

a 

0 

5 

12.210 

55.790 

143.917 

278.314 

455.985 

3 

a 

£ 

6 

11.208 

51.205 

131.761 

253.865 

419.438 

7 

10.411 

47.566 

122.267 

234.973 

387.374 

w 

1 

M 

\ 

8 

9.760 

44.592 

114.562 

219.847 

362.437 

1 

h 

B 

% 

9 

9.223 

42.108 

108.144 

207.421 

340.663 



10 

8.748 

39.989 

102.689 

198.878 

323.080 

Piecevise 
Continuous 
Model (N=3) 

8.776 

39.993 

101.758 

204.243 

329.268 


5. CONCLUDING REMARKS 

This paper proposed a piecewise continuous Timoshenko beam model which is to be used 
for the dynamic analysis of large flexible tapered beam-like structures. The procedure for 
establishing natural frequency has been described in detail. A tapered beam is divided into several 
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segments of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite 
element analysis,*the elosed-form solution of the Timoshenko beam equation has been used. 
Application of transfer matrix method relates all the elements as a whole. By corresponding 
boundary conditions and compatible conditions a characteristic equation for the global tapered 
beam has been yielded. Through the root-searching process to the characteristic equation the 
natural frequencies have been derived. A computer simulation is shown in this paper, and 
compared with the results obtained from the finite element analysis. While the comparable results 
is obtained, piecewise continuous Timoshenko beam model decreases the number of elements 
significantly. 
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